Surface Water-Groundwater Connections as Pathways for Inland Salinization of Coastal Aquifers.

Coastal agricultural zones are experiencing salinization due to accelerating rates of sea-level rise, causing reduction in crop yields and abandonment of farmland. Understanding mechanisms and drivers of this seawater intrusion (SWI) is key to mitigating its effects and predicting future vulnerability of groundwater resources to salinization. We implemented a monitoring network of pressure and specific conductivity (SC) sensors in wells and surface waters to target marsh-adjacent agricultural areas in greater Dover, Delaware. Recorded water levels and SC over a period of three years show that the mechanisms and timescales of SWI are controlled by local hydrology, geomorphology, and geology. Monitored wells did not indicate widespread salinization of deep groundwater in the surficial aquifer. However, monitored surface water bodies and shallow (<4 m deep) wells did show SC fluctuations due to tides and storm events, in one case leading to salinization of deeper (18 m deep) groundwater. Seasonal peaks in SC occurred during late summer months. Seasonal and interannual variation of SC was also influenced by relative sea level. The data collected in this study data highlight the mechanisms by which surface water-groundwater connections lead to salinization of aquifers inland, before SWI is detected in deeper groundwater nearer the coastline. Sharing of our data with stakeholders has led to the implementation of SWI mitigation efforts, illustrating the importance of strategic monitoring and stakeholder engagement to support coastal resilience.

Stormwater drives seasonal geochemical processes beneath an infiltration basin.

Deicing salt is an important component of road safety during winter storms. Stormwater infiltration best management practices aim to prevent the salt from polluting streams and waterways, but this may shift pollutants to groundwater resources. In response to limited field studies investigating groundwater quality impacts caused by input of salt from stormwater infiltration best management practices, we monitored water levels and quality of groundwater at various depths in an unconfined aquifer around a stormwater infiltration basin using in situ sensors coupled with grab sampling. Our observations revealed differences in groundwater chemistry with depth in the aquifer and processes that were driven by the seasonal changes in the chemistry of stormwater (salt-impacted in winter and fresh in non-winter) recharging the aquifer. Water-matrix interactions in the vadose zone beneath the basin affected the transport of sodium (Na) into groundwater following non-winter recharge. Sodium movement through the aquifer was delayed relative to chloride (Cl), indicating a longer residence time of Na in the vadose zone. Radium (Ra) concentrations were correlated with Cl concentrations, suggesting salt-impacted recharge caused desorption of Ra into groundwater because of increased salinity. Stormwater-influenced groundwater followed a preferential flow path due to heterogeneity of the aquifer materials, and water chemistry varied with time and location along the flow path. These results highlight the importance of well screen length, placement and depth, and frequency of observations when designing a monitoring network.

Hydrophysical and Hydrochemical Controls of Cyanobacterial Blooms in Coursey Pond, Delaware (USA).

Noxious cyanobacterial blooms are common in many ponds in the mid-Atlantic Coastal Plain. In Delaware, blooms normally occur between July and October, yet no in-depth analyses of the causes and predictors exist. A study using commercially available, high-frequency, continuous, and automated biogeochemical sensors at Coursey Pond, Delaware, a pond known for perennial summer blooms, was conducted to investigate how hydrophysical and hydrochemical conditions affect bloom dynamics. Cyanobacterial abundance (based on the in vivo phycocyanin fluorescence and phycocyanin/chlorophyll fluorescence ratios) increases during periods of high water temperatures (up to 32°C), low discharge through the pond (mean hydraulic residence time ≥5 d) with evaporative concentration of dissolved solids, and decreasing NO concentrations (reaching <0.1 mg L, the detection limit). These conditions lead to the uptake and depletion of bioavailable N in the pond surface waters and provide a competitive advantage for temperature-tolerant and N-fixing cyanobacteria. Irrigation water use within the watershed can exceed pond discharge during drier summer months, enhancing bloom-forming conditions. Bloom intensity varies due to storms but persists until mid-October to mid-November when temperatures cool, daylength decreases, and discharge and NO concentration recovers. Changes in these easily monitored physical and chemical parameters can serve to anticipate the onset, intensity, persistence, and the eventual dissipation of cyanobacterial blooms at Coursey Pond and similar ponds elsewhere. Therefore, the use of sensors and high-frequency data has the potential to assist in forecasting and management of blooms.

Assessing Coastal Plain Risk Indices for Subsurface Phosphorus Loss.

Phosphorus (P) Index evaluations are critical to advancing nutrient management planning in the United States. However, most assessments until now have focused on the risks of P losses in surface runoff. In artificially drained agroecosystems of the Atlantic Coastal Plain, subsurface flow is the predominant mode of P transport, but its representation in most P Indices is often inadequate. We explored methods to evaluate the subsurface P risk routines of five P Indices from Delaware, Maryland (two), Virginia, and North Carolina using available water quality and soils datasets. Relationships between subsurface P risk scores and published dissolved P loads in leachate (Delaware, Maryland, and North Carolina) and ditch drainage (Maryland) were directionally correct and often statistically significant, yet the brevity of the observation periods (weeks to several years) and the limited number of sampling locations precluded a more robust assessment of each P Index. Given the paucity of measured P loss data, we then showed that soil water extractable P concentrations at depths corresponding with the seasonal high water table (WEP) could serve as a realistic proxy for subsurface P losses in ditch drainage. The associations between WEP and subsurface P risk ratings reasonably mirrored those obtained with sparser water quality data. As such, WEP is seen as a valuable metric that offers interim insight into the directionality of subsurface P risk scores when water quality data are inaccessible. In the long term, improved monitoring and modeling of subsurface P losses clearly should enhance the rigor of future P Index appraisals.

Model evaluation of denitrification under rapid infiltration basin systems.

Rapid Infiltration Basin Systems (RIBS) are used for disposing reclaimed wastewater into soil to achieve additional treatment before it recharges groundwater. Effluent from most new sequenced batch reactor wastewater treatment plants is completely nitrified, and denitrification (DNF) is the main reaction for N removal. To characterize effects of complex surface and subsurface flow patterns caused by non-uniform flooding on DNF, a coupled overland flow-vadose zone model is implemented in the multiphase flow and reactive transport simulator TOUGHREACT. DNF is simulated in two representative soils varying the application cycle, hydraulic loading rate, wastewater quality, water table depth, and subsurface heterogeneity. Simulations using the conventional specified flux boundary condition under-predict DNF by as much as 450% in sand and 230% in loamy sand compared to predictions from the coupled overland flow-vadose zone model, indicating that simulating coupled flow is critical for predicting DNF in cases where hydraulic loading rates are not sufficient to spread the wastewater over the whole basin. Smaller ratios of wetting to drying time and larger hydraulic loading rates result in greater water saturations, more anoxic conditions, and faster water transport in the vadose zone, leading to greater DNF. These results in combination with those from different water table depths explain why reported DNF varied with soil type and water table depth in previous field investigations. Across all simulations, cumulative percent DNF varies between 2 and 49%, indicating that NO3 removal in RIBS may vary widely depending on operational procedures and subsurface conditions. These modeling results improve understanding of DNF in RIBS and suggest operational procedures that may improve NO3 removal.